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The Density Contrast of the Shapley Supercluster 
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ABSTRACT 

We calculate the density contrast of the Shapley Supercluster (SSC) based on 
the enhanced abundance of X-ray clusters in it using the extended Press-Schechter 
formalism. We derive a total SSC mass of Mtot = (4.4 ± 0.44) x 10^^ Mq within a 
sphere of 50Mpc centered at a distance of about 160 Mpc. The nonlinear fractional 
density contrast of the sphere is (1 + (5) = 1.76 ± 0.17 relative to the mean matter 
density in the Universe, but the contrast increases in the interior of the SSC. Including 
the cosmological constant, the SSC region is found to be gravitationally unbound. The 
SSC contributes only a minor portion (9.0% ± 2.1%) of the peculiar velocity of the 
local group. 
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1 INTRODUCTION 



The Shapley Supercluster (SSC) is one of the l argest 
known structures in the local Unive r se l[Fabianl flQQll: 
Quintana et al. 19951: Ettori et al.lll997 : Reiseneeeer et al.l 



2000l: 



de Filippis et al 
Proust et al.1 120061 ). Recent 



20051 : iHaines et alT l2006l : 

measurements based on 
the TuUy-Fisher relation give 11645 ± 859 km s^^ for the 
Hubble flow recession velocity of A3558, one of the central 
clusters in the SSC (|Springob et al.l l2007h . This corre- 
sponds to a redshift of z = 0.0388, a luminosity distance 
of (166 ± 13) Mpc, and an angular diameter distance of 
(154 ± 11) Mpc. Despite the great distance of the SSC, 
it has recently been argued that the SSC contributes 
significantly to the peculiar velocity of the local group (LG) 
based on the enhanced ab undance of X-ray clusters in it 
jKocevski fc EbeUndlioO^ ). 

The contribution of the SSC to the LG peculiar veloc- 
ity depends critically on the mean matter density in the su- 
percluster. It is a non-trivial matter to relate the enhanced 
abundance of massive X-ray clusters in the region to its over- 
density, since clusters are strongly bia sed with respect to 
the underlying d i stribution of matte r (IMo fc White 1 19961 : 
ISheth et all l2(30ll : lEvrard et aP I2OO2I : iBahcall et al.l l2004h . 
Because X-ray clusters correspond to rare density peaks, 
even mild enhancements in the mean matter density can 
result in the formation of many more clusters within the su- 
percluster. Since this bias depends on the cluster mass, the 
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relationship between a sample containing objects in a range 
of masses and the underlying matter cannot be described 
accurately by a constant bias factor. 

In this paper we study the dependence of the num- 
ber of X-ray clusters in the SSC on the underlying matter 
overdens ity in this region . We u se the hybrid model pro- 
posed bv lBarkana fc Loebl l|2004h - combining the ext ended 
Press-Schechter (ePS) formalism (iBond et al. 1991) with 
the st ructure formation prescription of Isheth fc TormenI 
(|l999l ) (ST). The model describes well published results 
from numerical simulations (see discussion around Fig. 4 of 
Barkana & Loeb (2004)). Using the mass-dependent, non- 
linear bias for X-ray clusters, we calculate the matter over- 
density in the SSC region from the number of observed X-ray 
clusters there. 

In 921 we describe how the iBarkana fc Loebl (|2004| ) hy- 
brid model can be used to calculate the mass in a large 
region enclosing collapsed objects. We then specifically con- 
sider the SSC in [JSJand describe the sample of x-ray clusters 
used in our calculation. In ^we apply our formalism and 
calculate the matter overdensity, cluster overdensity, cluster 
bias, and total mass in the SSC region. Given these results, 
we then use the spherical collapse model in f|5] to consider 
the dynamics of the region. In fj6l we estimate the SSC con- 
tribution to the peculiar velocity of the LG. In SJT] we explore 
the radial dependence of our results in the SSC and, finally, 
investigate the robustness of our results in |8] 



We assume 
ogy with the 
cal parameters 



a flat, ACDM cosmob 
standard set of cosmologi- 

{Q,M.o,^A,o,^b.o,h,(Ts.,ct,r) = 
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(0.279,0.721,0.0462,0 .701, 0.817, 0.960, 0.000) 
jKomatsu et al.ll2008l ). 



2 METHOD 



According to the ePS prescription l|Bond et al.lll99ll ). if the 
linear density fluctuations in the universe are extrapolated 
to their values today and smoothed on a comoving radius R, 
a point whose overdensity exceeds a critical value of Sc{z) ~ 
1.68 D{0) / D{z), belongs to a collapsed object with a mass 
M — (4/3) IT pcritfi ^M,o if R is the largest scale for which 
the criterion is met. Here, D{z) is the linear growth factor at 
redshift z, pcrit,o is the closure density of the Universe, and 
f^Ai.o is the matter density parameter today. The critical 
value of the overdensity is also known as the barrier. For 
a Gaussian random field of initial density perturbations, as 
indicate d by measurements of cosmic microwave background 
(CMB) (|Komatsu et ahlbOOSl ), the probability distribution 
of the extrapolated and smoothed over-density, 5r, is also a 
Gaussian: 



dP{SR,SiR)) 
dSn 



dSi, 



V^T^SiR) 



exp 



2S{R) J 



with zero mean and a variance given by: 



S{R) 



(1) 



(2) 



where P (fc) is the linear power-spectrum of density fluc- 
tuations today as a function of comoving wave-number k, 
and kmax = 1/R. Since equation Q is a monotonically de- 
creasing function of R (or M), the smoothing scale can be 
uniquely specifled by the variance of the over-density field 
smoothed on that scale. 

If we define f{Sc{z),S)dS to be the fraction of mass 
contained in halos in the mass range [M, M + dAI] corre- 
sponding to [5*, S + dS] at redshift z, then the comoving 
number density of collapsed objects in that mass range is 



dn 
dM 



f^M.O pcritfi 
M 



dS 



dM 



f{Sc{z),S). 



(3) 



The ePS prescription, which uses a constant barrier indepen- 
dent of mass with a value, Sc{z), derived from the spherical 
collapse model, gives 



fpsi5ciz),S) dS 



/27 S 



exp 



dS, 



(4) 



where u = 5c{z)/y/S is the number of standard deviations a 
density fluctuation on a scale S must be above the mean to 
have crossed the critical density threshold. Incorporating a 
moving (i.e. scale dependent) barrier generated from ellip- 
tical c ollapse with two free p arameters, the ST prescription 
gives l|Sheth fc Tormenlll999l ) 



/st(5c(z),S) dS = A 



S 727 



(a V 



2\p 



exp 



a V 



dS 



(5) 

and results in a mass function that better matches simu- 
lations. The best-fit values of the parameters are a = 0.75 
and p = 0.3, while the normalization factor is y4 = 0.322 
ijBarkana fc Loebll2004l '). 

The unconditional mass function generated with 



f{Sc{z),S) represents an average over all regions of space 
(or equivalently, over all realizations of density Fourier mode 
amplitudes). It assumes no prior knowledge of the overden- 
sity on a given scale. If we fix the average linear overdensity 
to be 5l in a region smoothed on a particular scale. So, we 
can generate a conditional mass function with a conditional 
form of / = f{5c{z),S\5L,So). In ePS, this can be done by 
substituting 5c{z) — > 5c{z) — Sl and S — > S — So and gives 
results that agree well with simulations. The same substitu- 
tion with fsT, however, does not agree as wel l. This is be- 
cause the new barrier height is scale dependent (|Zhang et al.l 

:2ood)- 

However. iBarkana fc Loebl l|2004h suggested using a hy- 
brid model in which the conditional mass function is gen- 
erated from contributions by both fps and fsT in regimes 
where they each fit best with simulations. The resulting col- 
lapsed mass fraction per variance interval is 



tp.s{Oc{z),S) 

(6) 

This model describes well the numerical resu lts from cos- 
mological simulations l|Barkana fc Loebl |2004| ). Because < 
5l >~ in any region, the average number of collapsed 
objects, iV', with mass corresponding to 5* in a region 
corresponding to So is generated by equation @ with 
fhybrid{Sc{z), S\5l = 0, So). We use this prescription to map 
values of the overdensity in a region to the average number 
of collapsed objects contained in it. Variations between this 
conditional average number of objects, N and the actually 
number, A'', of those counted in the region result only from 
Poisson fluctuations. The cosmic variance has been taken 
out by stipulating S = Sl in S = Sq. 

Given an integer number, A^, of observed objects resid- 
ing in a region So, the differential probability distribution 
of values of N that have resulted in N is. 



dP{N\N) 
dA 



= B Ppor. 



dPjSL) dSL 

dSL dN ' 



where 



Ppoisson{k, A) 



fc! 



(7) 



(8) 



is the Poisson distribution, dP{5L) / dSL is the unconditional 
distribution of overdensities in the region So given by equa- 
tion m. The Jacobian dSL/dN is derived from equations ^ 
and (ini, and the coefficient B is set so as to normalize the 
integral of equation (O over A' to unity. 

When applying this formalism to observations, there is 
an additional complication in that overdensities and the sizes 
of regions are observed in Eulerian coordinates which evolve 
as the region breaks away from the Hubble fiow, while the 
ePS and ST prescriptions rely on initial values in Lagrangian 
coordinates. Lagrangian sizes in the comoving frame do not 
change over time. As in Im unoz fc Loebl (|2008l ). we use the 
spherical evolution model to calculate the Lagrangian size, 
Rl, corresponding to the observed Eulerian size. Re, of a 
region containing the linear overdensity, 5l in a KCDM uni- 
verse. The extent of the collapse depends on the magnitude 
of the overdensity. The more overdense the region is, the 
larger it would have to be initially in order for it to col- 
lapse to the same value of Re. Similarly, a lower value of 
the overdensity would mean that the material inside Re 
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Table 1. X-ray Cluster Sample 



Cluster 


RA;Dec 
(J 2000.0} 


z 


r 


fx 

lA n— 12 

(10 erg 


k T 


Ml625 


Mhalo 


Refer. 




(deg) 




(Mpc) 




{keVj 


/I nl4 n ^ \ 
(10 M0) 


/I nl4 Ti/T \ 

(10 Mq) 




B6 


1 Q4 YQi^- 91 Q1 1 




33 


3.66 


1.8 


89 


1.9 




A3548 


1 Q8 37() -44 Qfi3 




39 


4.26 


1.9 


0.97 


2.0 


A 


Bl 


1 Qfi oj^n- 1 7 nni 




44 


7.53 


2.5 


1.5 


3.0 


ji^ 


A721S 


1 Qfi f^il '\7 fizL9 


n n4Qn 

u. uiyu 


23 


7.76 


2.5 


1.5 


3.0 


ji^ 


PT'ZA 11410 zL494fi 


91 9 R1 Q- 49 777 


u. uiyu 


41 


9.35 




1.7 


3.5 


Q 


A1631 


1 Q'\ 949- 1 fi '^70 








9 s 
z.o 


1 7 






A1736 


901 7^R- 97 1 '^'^ 


0.0453 


17 


27.54 


3.0 


1.9 


4.0 






9n'i 1 HQ- 9<A 9 




16 


11.90 


3.0 


1.9 


4.0 




3528S 


193.673:-29.231 


0.0528 


31 


12.20 


3.1 


2.0 


4.2 


A 


A3530 


193.917:-30.367 


0.0537 


33 


9.24 


3.2 


2.1 


4.4 


A 


A3528N 


193.598:-29.010 


0.0528 


31 


10.53 


3.4 


2.3 


4.8 


A 


RX J1252.5-3116 


193.143:-31.266 


0.0535 


33 


16.09 


3.8 


2.7 


5.7 


A 


SC 1327-312 


202.514:-31.664 


0.0495 


6.6 


12.25 


3.8 


2.7 


5.7 


A 


A3556 


201.001:-31.656 


0.0479 


2.5 


1.72 


3.8 


2.7 


5.7 


A 


A3528 


193.640:-29.129 


0.0528 


31 


24.32 


4.0 


3.0 


6.2 


A 


SC 1329-314 


202.875:-31.812 


0.0446 


15 


5.84 


4.2 


3.2 


6.7 


A 


A3562 


203.446:-31.687 


0.0490 


5.6 


29.16 


4.3 


3.3 


6.9 


B 


A3532 


194.336:-30.375 


0.0554 


38 


21.35 


4.4 


3.4 


7.1 


A 


A1644 


194.332:-17.381 


0.0473 


45 


4.45 


4.6 


3.7 


7.6 


B 


A3558 


202.011:-31.493 


0.0480 


0.0 


57.87 


4.9 


4.0 


8.4 


B 


A3571 


206.860:-32.850 


0.0391 


40 


110.9 


6.8 


6.6 


13.7 


B 



Col. (1): Source name. Col. (2): RA and Dec (J2000). Col. (3): Redshift. Col. (4): Distance from A3558 in Mpc. Col. (5); /x in the 
0.1 — 2.4 keV energy band in units of lO"^'^ ergs cm~-^ s~^ . Col. (6): Cluster average temperature in keV. Col. (7): M1625 calculated 
from equation HIOI I in units of Mq. Col. (8): Mf^aio = 2.1Mig25 in units of Mq. We conservatively estim ate errors in Mi^gio a nd Mig25 
to be about 3 0%. Col. (9): Literature reference for the cluster temperature: (A) Ide Filippis et al.l ||2005| '): (B) IVikhliiiir] 1I2OO8I '): and (C) 
lEbeling et al] f2002i) . 

function given by ePS and similar prescriptions ignores sub- 
structure, we must be careful to consider only the largest 
structures. Fortunately, the largest collapsed halos are also 
the most sensitive tracers of the mass function, being on its 
exponential tail. We consider the most massive. X-ray lumi- 
nous clusters in the supercluster as tracers of the most mas- 
sive halos and make the assumption that each halo hosts one 
such cluster. This is a good approximation given how well 
the average ST halo mass function in the universe matches 
the cluster mass function (jVikfilinin 2 00S ). Our tendency to 
focus on only the rarest objects is moderated by the need to 
reduce Poisson fluctuations on our sample. 

To balance these requirements, we construct a sam- 
ple of clusters in the SSC whose host halos have masses 
above Mmin ~ 1.75 x 10^* Mq. The sample consists pri- 
marily of clusters studied by Ide Filippis et al.l (2005), but 
contains an additional cluster from the Clusters i n the Zone 
of Avoidance (CIZA) sample (|Ebehng et al.ll2002l ') . We have 
approximated the distance between each of the clusters in 
our sample and the nominal center of the SSC, which we 
have taken to be A3558. To estimate these distances, we 
assume that each clu ster has the same pec uliar velocity as 
A3558, measured by ISpringob et af] (|2007l ). This is nearly 
equivalent to assuming no relative peculiar velocity among 
the clusters. As we will see in !}S]this approximation is ad- 
equate for our purposes. Our entire sample extends to a 
radius of roughly Re ~ 50 Mpc around A3558. At a dis- 
tance of 164 Mpc, 50 Mpc perpendicular to the line-of-sight 
corresponds to ~ 1 7.5 degrees. We a s sume that the clusters 
newly detected by Ide Filippis et al.l l|2005h have z = 0.048. 



came from a relatively smaller Lagrangian size. Integrating 
the equations of motion results in a mapping between the 
comoving Eulerian size of a viewed region and the comoving 
Lagrangian size of the region from where the same material 
originated in the early universe. For a fixed value of Re, 
there is a one-to-one relationship between Rl and the value 
of 5l that collapses Rl to Re. Rl is then associated with 
5-0. 

Since matter shells in the spherical collapse model do 
not cross until collapse, the amount of matter inside Rl is 
the same as that inside Re. The mass contained within the 
observed size Re is 

Mtot{< Re) = i-i/S) TV rIQm.o Pcr,t,o {1 + Si), (9) 

where the initial value of the overdensity is 5i = 
D{zi) 5l/D{zo)), Zi is the initial redshift before the re- 
gion begins to evolve nonlinearly, and zq is the observed 
redshift of the region. The nonlinear matter overdensity 
is then S — M/{Vssc ^m,o Pcrit.o) — 1, and the nonlinear 
bias is & = {N/N' ~ 1)/S, where Vssc = (4/3)7ri?| and 
N' — nsrVssc (Rl/Re)^ . nsT is the unconditional, ST 
mass function. 



3 THE CLUSTER SAMPLE 

To calculate the mass and overdensity described in the pre- 
vious section for the SSC, we need to know the number of 
the collapsed objects in the SSC, their minimum mass, and 
the size of the region in which they reside. Since the mass 
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We make the same assu mption about ^3548. We exclude the 
object denoted as B8 bv lde Fihppis et all l|2005h because it 
is not a "confirmed cluster". Table[l]gives our resulting sam- 
ple of 21 clusters. It includes a list of the clusters, positions, 
redshifts, estimated distances to A3558, measurements of 
the X-ray flux and temperature {fx and Tx, respectively), 
and the resulting estimates of the host halo mass, M1625 
and Mhaio, which are define d below. Our sample is similar 
to the 17 clusters studied bv lKocevski fc Ebelind |200Q) in 
the SSC in number and extent. 

Here, M1625 is the mass of the halo out to the radius 
within which the matter density is 1625 times the mean 
matter densit^Q. To calculate M1625 , we make use of its re- 
lationship to X-ray temperature Tx (|Vikhlininll2008h : 



Mi625 = (2.95 ±0.10) X 10^ 



fca Tx 
5keV 



(10) 

For CIZA J1410.4-42 46, we calculate M1625 from the X-ray 
flux via the relation (|Vikhlininll2008l ): 

InLx = (47.392 ± 0.085) -f (1.61 ±0.14) ZnMi625 
±(1.850 ± 0.42) In E{z) ~ 0.39 ln{h/0.72) 
±(0.396 ±0.039), (11) 

where the X-ray luminosity, Lx — {iii dj^) fx , is in units of 
ergss"^, (11 is the luminosity distance, M1625 is in units of 
Mq, and E{z) = y/{nM,o {1 + z)^ + ^a). 

Mhaio, on the other hand, is the mass for which ePS-like 
prescriptions generate the mass function and corresponds 
roughly to the mass out to the radius within which the mat- 
ter density is 180 times the mean density. The relationship 
between various defin itions of mass can easily be calculated 
||Hu fc Kravtsovll2003l ) . Values of M1625 underestimate those 
of Mhaio by ~ 52% (i.e. Mhaio = 2.1Mi625)- We estimate 
errors in these masses to be about 30%, resulting primarily 
from error in the temperature measurement and intrinsic 
scatter in the AI — Tx and M — Lx relations from equa- 
tions ^ and (fTTI) . 

Because the clusters under consideration are very mas- 
sive and luminous, we ignore the effect of the cluster se- 
lection function. The flux of our dimmest cluster is bright 
enough to be seen even in regions where ROSAT sensitiv- 
ity is slightly lower (i.e. in the ~ 1.5 deg patch centered at 
(RA:DEC) = (199.295:-34.393)). However, absorption by the 
galaxy could hide a few additional clusters. On the other 
hand, because the core of the SSC is clearly more overdense 
than the region out to 50Mpc, our assumption that each 
cluster represents one halo may not hold there. Multiple 
clusters in the core may lie within a single halo (or multiple 
halos in the process of merging). Thus, the number of ha- 
los may be somewhat different than the 21 we assume. Our 
analysis, however, is independent of the individual masses 
of these halos and depends only on the minimum halo mass 
of our sample and the total number of halos. Because the 



^ The notation M500 used sometimes in the literature to refer to 
the halo mass out to the radius containing an a verage matter de n- 
sity equal to 500 times the critical density (eg. I Vikhlininl |200^). 
For our cosmology and at the redshift of the SSC, this notation 
corr esponds to M16 2 5 use d here and elsewhere in the literature 
(eg. iHu fc Kravt"sOTl | |2003D ') where the average matter density in 
the sphere is 1625 times the mean matter density of the universe. 



error in our mass measurements is only 30 %, errors in the 
minimum mass can be described by additional error in the 
number of clusters being above a fixed threshold. In addi- 
tion, we expect the variation in A*' due to errors in Mhaio to 
be small. If it turns out that all of our mass estimates are 
systematically 30 % too high, then only two clusters below 
the limit would have made it into our sample undeserved. 
In !j8l we show the dependence of our results on the exact 
value of N. 



4 SHAPLEY OVERDENSITY 

In Figure [T] we plot the distributions of several quantities 
calculated through the prescription outlined in i|2]or derived 
from the resulting quantities. Panels (a)-(i) show results for 
the expected number of collapsed objects above Mmin, the 
linear matter overdensity, 5l, the total mass of the SSC, 
Mtot, the nonlinear matter overdensity, 5, the cluster over- 
density, Sn, the cluster bias, b, the Lagrangian radius of the 
SSC, Rl, the contribution of the SSC to the peculiar ve- 
locity of the LG, Vp, and the peculiar radial velocity of the 
outer edge of the SSC, Vrad,p- 

The distribution of N in panel (a) shows a mean clearly 
less than our assumed value of N — 21. This is due to the 
[dP/dSL) (dSL/dN) factor in equation ([T]). If, in the absence 
of information about N, each value of A^^ were equally likely, 
we would have < A" >= N. However, including the true 
prior probability distribution means that lower values of N 
are more likely. 

The linear overdensity can be interpreted given the 
value of cr = on the scale of the supercluster. For 
R = 50Mpc (64Mpc), a « 0.23 (0.18). The mean lin- 
ear overdensity of ~ 0.47 in panel (b) represents a ~ 2.0a 
(~ 2.6a) fluctuation in the overdensity in the region. As a 
sanity check, we note that the probability of a linear den- 
sity as high or higher, assuming initial Gaussian fluctua- 
tions, is approximately 1/50 (1/200), and there are about 
[200Mpc/50Mpc]3 ~ 60 ([200 Mpc/60 Mpc]^ ~ 40) such re- 
gions out to the distance of the SSC. So it is not unlikely 
that we do observe such a region within 200Mpc, nor is it 
so unlikely that the supercluster remains one of the most 
overdense regions out to that distance. 

Our analysis gives Mtot = (4.4 ± 0.44) x 10^® for 
the mass of the SS C. This is similar to th e mass in the 
region estimated by Ide Filippis et al.1 (|2005|) and o nly a 
bit lower than esti mates by Reisenegger et ah! l|2000l) and 
iProust et all (|2006i). 

Panels (d)-(f) show the probability distributions of the 
matter overdensity, S, the cluster overdensity, 5n, and the 
cluster bias, b. The nonlinear overdensity distribution shows 

5 = 0.76 ± 0.17 and Sn = 2.4 ± 0.32. Larger values of S 
correspond to lower values of Sn and b since the number of 
clusters is fixed and so more matter in the region gives closer 
agreement between the matter and cluster densities. The 
bias is 6 = 3.5±1.4, however, the distribution is not Gaussian 
and is skewed toward values higher than its peak at 6 « 1.6. 
This peak value is just over 1 — a hi gher than the maxi- 
mum of the range of 6 = 1.2-1.9 that iKocevski fc Ebelinj 
(|2006l ) estimated (for Qm,o = 0.279) by assuming that their 
calculation of the peculiar velocity of the LG induced by 
their sample is equal to the true value. We calculate only 
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Figure 1. Probability distributions of the expected number of collapsed objects above Ai^nin, the linear matter overdensity, 5^, the total 
mass of the SSC, Mtot, the nonlinear matter overdensity, <5, the cluster overdensity, 5pf, the cluster bias, b, the Lagrangian radius of the 
SSC, Rl, the contribution of the SSC to the peculiar velocity of the LG, Vp, and the peculiar radial velocity of the outer edge of the 
SSC, v^ad,p- E'nd Vrad,p plotted in units of kms~^. The distributions assume N = 21, M^in = 1-75 X lO-*^* A/q, and = 51 Mpc. 



a 6% chance that the bias falls into the b = 1.2-1.9 range. 
For reference, we find that b = 1.5 corresponds to 5 ~ 1.2. 
The discrepancy in the bias may result from differences be- 
tween the clusters a s sociat ed with the SCC and the entire 
iKocevski fc Ebelind (|2006l ) sample. The entire sample in- 
cludes objects that are less intrinsically luminous (i.e. less 
massive) than those in the SSC. 



5 SHAPLEY SUPERCLUSTER DYNAMICS 

While we used the spherical collapse model in previous sec- 
tions to calculate Rl corresponding to Re for each value of 
5l , we now look at the results from the model in detail about 
the dynamics of the SSC region. We describe our analysis in 
ij5.ll and possible tests for our results in H5.2I 

5.1 The Spherical Collapse Model 

The evolution of the physical size r of a region containing 
mass Mtot is described by: 

_ /_l_y _ / 1 1 dH\ ds 

where Mtot is the mass in the region given by equation @, 
s — r/vi, r is the evolving Eulerian size of the region, and 
where the evolution of the Hubble parameter is given by the 



Friedmann equation, 

H^^^=HoJ^+nA.o. (13) 
a at \ a-' 

The evolution of the matter and dark energy density param- 
eters {Qm and SIa, respectively) are given by 

nM = (i-nA)^ „ ^3 . (14) 

Subscripts "0" refer to values today, and subscripts "i" re- 
fer to initial values at at. We take ai — 0.01, set Si = 1 
(i.e. r — ri) and {ds/da)i = (1 — (5i/3) (l/tsi) as our initial 
conditions, and numerically integrate equation (|12p up to 
a/ = 1/(1.0388) ^ 0.962, corresponding to the Hubble flow 
redshift of A3558. 

In setting (ds/da); = (1 — 5i/3) (l/a^), as appropriate 
for a growing-mode perturbation, we have assumed that the 
region is not initia lly expanding with the Hubble flow as in 
iLokas fc Hoffman! ([2001). The —6i/3 term accounts for the 
peculiar velocity of the region due to its overdensity. Since 
we expect the Si to be small, however, the peculiar velocity 
term will also be small. 

The probability distribution for the Lagrangian radius, 
Rl = ri/ ai, is given in panel (g) of Figure [T] The distribu- 
tion gives Rl = (64 ± 2.0) Mpc, indicating that the region 
of the SSC deviates from the Hubble flow and collapses by 
only ~ 20%. 

The more interesting result from applying the spheri- 
cal collapse model to the SSC is that the region is still ex- 
panding at a radius of 50 Mpc. While not expanding as fast 
as the Hubble flow, the SSC has not reached turn-around. 
Moreover, we find that, due to the repulsive effect of the cos- 
mological constant, the region will never reach turn-around 
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Figure 2. The evolution of spheres around the center of the SSC 
as a function of the scale factor, a, calculated from the spheri- 
cal collapse model. Results are shown for radii of 51 Mpc (solid 
line), 20 Mpc (dot-dashed line), 10.3 Mpc (long-dashed line), and 
10 Mpc (short-dashed line). Since a = 1 today, values of a > 1 
refer to time in the future. The vertical line denotes the redshift 
of the SSC. The mean values of the overdensity within the radii 
shown have been adopted in plotting the trajectories. 



despite its overdensity but will continue to expand forever. 
This can be seen in the trajectory of the 51 Mpc region plot- 
ted in Figure [2] Trajectories have also been plotted for the 
interior regions of the SSC (see SJT]). 

We arrive at the peculiar radial velocity of the region 
by subtracting the Hubble velocity from the total expansion 
rate. The distribution in the peculiar radial velocity of the 
50 Mpc shell, Vrad,p, is given in panel (h) of Figure [T] and 
gives Vrad.p = (— 437± 87) kms"'^. The minus sign indicates 
that the region is expanding somewhat slower than the ~ 
3600 km s~^ Hubble velocity at that radius. 



5.2 Observational Tests of the Model Predictions 

Our prediction for the outward velocity of the 50 Mpc shell 
could potentially be tes ted with the TuUy-F i sher relation or 
by T y pe la supernova e (|Masters et al.ll2006l : ISpringob et al.l 
I2OO7I : IJha et al.|[2007l ). By determining the distance to an 
object of a known redshift, one may infer its peculiar ve- 
locity by subtracting oflF the Hubble velocity at its inferred 
distance. The quantity Vrad,p is the radial peculiar velocity 
with respect to the center of the SSC. Thus, the peculiar ve- 
locity of the SSC itself must be accounted for when using the 
Tully -Fisher relation to m easure v,.ad,p- The latest results 
from ISprineob et"all (|2007l ) provide a peculiar velocity for 
A3558, the central cluster of the SSC, of (2980±859) kms'^ 
However, the error on the peculiar velocity of A3558 is too 
large to constrain Vrad,p reliably. Below we explore alterna- 
tive observational methods for future measurements of the 
predicted peculiar velocities. 



The peculiar velocity of an X-ray cluster can be directly 
inferred from its contribution to the CMB through the ki- 
netic Sunyaev-Zel'dovich (kSZ) effect in which CMB pho- 
tons scatter off free elections moving along with the bulk 
velocity of the cluster. The temperature fluctuation induced 
in this way by a cluster is given by: 



AT V 



(15) 



where t ~ axrie Rc is the optical depth for Thomson scat- 
tering of the CMB by free electrons in the cluster, ctt is the 
Thomson cross-section, rie is the number density of elec- 
trons in the cluster, Rc is its core radius, and v is its pecu- 
liar velocity. The minus sign indicates that velocities away 
from us result in temperature decrements. Because A1644 
is on the edge of the SSC (albeit the edge on the sky as 
opposed to along the line of sight), we use it as a test case 
to estimate AT/T due to Vrad.p of a hypothetical cluster 
on the 50 Mpc shell and positioned along the line of sight 
to A3558. The radius of A1644 is ~ 0.7 Mpc, while its X- 
ray luminosity is L x ~ 10^"^ ergs s"'^ and its temperature 
is about 6 x lO^K (|de Filippis et al.l [20051 ). Accounting for 
the X-ray luminosity through thermal bremsstrahlung we 
get ne ~ 5 X 10~''cm~'^. If the hypothetical cluster under 
consideration (whose properties are the same as A1644) is 
50 Mpc in front of A3558 directly along the line of sight so 
that Vrad,p is also along the line of sight, the CMB tempera- 
ture fluctuation due to Vrad,p would be AT/T ~ —2 x 10~® 
with AT ~ —bfiK. Because the peculiar velocity is radially 
inward toward the center of the SSC, if the cluster is be- 
tween us and A3558, then this velocity is away from us and 
results in a temperature decrement. 

The SSC as a whole also induces CMB signals through 
the kSZ effect from to its own peculiar velocity and through 
the Rees-Sciama effect (Rccs & Sciama 1968i) . For an over- 
density of (5 = 0.76, the optical depth to electron scat- 
tering through the SSC is r = 2.3 x 10"^. The resulting 
kSZ signal due to the peculiar velocity of the structure is 
AT/T « -2.1 X 10"^. In the Rees-Sciama effect, the evolu- 
tion of the gravitational potential, </> = —GMtot/R, during 
the photon crossing time of the system introduces additional 
temperature fluctuations of order ~ If the poten- 
tial becomes deeper (i.e. more negative) as a photon passes 
through the region, the redshift experienced by the pho- 
ton on its way out of the potential well is greater than the 
blueshift it receives on its way in and the CMB in that di- 
rection shows a temperature decrement. However, between 
z = 0.051 and z = 0.027 (which spans the light crossing 
time for the SSC of ~ [lOOMpc/c]), the SSC grows from 
50.47Mpc to 51.54Mpc. For Mtot = 4:A x IO^H/q, the 
change in the potential resuhs in AT/T ^ 8.7 x 10~^. The 
kSZ and Rees-Sciama effects due to the SSC as a whole act 
against each other. 

The predicted signal can probed by a 
number of upcoming CMB observatories such 
as Planck (http://planck.esa.intl, the Ar- 
cminute Cosmology Bolometer Array Reciever 
( http: / / cosmology.berkeley.edu / group / swlh/acbar ) , 
the Atacama Pathfinder Experiment 

(http://www.mpifr-bonn.mpg.de/div/mm/apex I, the 
South Pole Telescope (http://spt.uchicago.edu|), 
and the Atacama Cosmology Telescope (ACT, 
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Figure 3. The radial profile of the SSC. The number of observed 
clusters, N , the matter overdensity, 5, and the amount of excess 
mass, A Mtot within as functions of are plotted in the top 
three panels. The bottom two panels show the peculiar velocity of 
the local group due to the overdensity contained within Re and 
the radial peculiar velocity of the shell with radius Re centered on 
A3558. Solid lines indicate mean values, while the dashed curves 
denote I — a bounds due to Poisson fluctuations. 



I ht t p : / / www .physics.princet on.edu/ act ) . With a tem- 
perature noise of ~ 35 K per pixel, WMAP is not 
sufficiently sensitive to detect this effect. On the other 
hand, ACT has a temperature noise of ST ~ 5 nK per 
pixel and a pixel size of 3 arcminutes^ . At a distance of 
(150 - 50 = 100) Mpc, an A1644-like cluster would take 
up ~ 600 pixels, lowering the detectable level to ~ 0.2 /iK. 
Therefore, the kSZ signal due to this cluster is, in principle, 
detectable with ACT at the ~ 26 a level. Though additional 
sources of error ma y make this detection more challenging 
(jSehgal et all '2005). the major obstacle to measuring this 
signal is disentangling it from the primary anisotropy that 
is at least an order-of-magnitude larger and has an identical 
spectrum. It is also important to note, that some type 
of distance measure from something like the TuUy-Fisher 
relation or from Type la supernovae is still necessary to 
determine which clusters are on the 50 Mpc shell. 



6 CONTRIBUTION TO PECULIAR 
VELOCITY OF THE LOCAL GROUP 

Next we calaculate the contribution of the Shapley Supler- 
cluster to the peculiar velocity of the LG, Vp, given our result 
for the SSC mass. The induced velocity in the linear regim43 



^ The expression for the peculiar velocity given in equation l|16|l 
is only valid in the linear regime. However, since the SSC over- 
density is modest, the formula gives a good approximation to the 
true value. 



is given by (iPeeblesl l 19931) 



47r 



(16) 



where / — (a/D) (dD/da) ~ and y and x point from 
the center of the SSC to each SSC mass element and to the 
LG, respectively. Because we do not explore the mass dis- 
tribution within the SSC, we assume a constant value of the 
overdensity out to 50 Mpc. The distribution of values for Vp, 
reflecting the distribution in the SSC overdensity, is plotted 
in panel (i) of Figure [1] We find Vp — (55 ± 13) kms~^, only 
(9.0 ± 2.1)% of the 612km s~^ peculiar velocity o f the LG 
with respect to the CMB l|Loeb fc NaravanI 120071). This is 
much less than the value estimated by Kocevski fc Ebelinej 
(200^) of ~ 30%. Even h = 1.5 and the corresponding value 
of 5 = 1.2 (see Q produce only Vp ~ 92kms~^, but this 
is almost 3 — a from our mean va l ue. Thus, we agree quali- 
t ative ly with Ravchaudhurv et al. l|]"99ll ). lReisenegger et al.l 
( 20001) ■ and lLoeb fc NaravanI (|2007l ) that the SSC does not 
contribute very significantly to the peculiar velocity of the 
LG. 



7 RADIAL PROFILE 

The number of clusters in our sample decreases sharply with 
radius from the center of the SSC, and the clusters are not 
numerous enough to adequately sample the entire volume 
of the SSC. Nonetheless, we can obtain some interesting es- 
timates about the radial distribution of the quantities we 
investigate by considering subsamples of clusters from Ta- 
ble[T]with distances from the center of the SSC within vary- 
ing values oi Re- 

Figure [3] shows the radial profiles we calculate for the 
number of observed clusters, the matter overdensity and the 
amount of excess mass within Re as functions of Re- Also 
shown are the peculiar velocity of the LG due to the over- 
density contained within Re and the radial peculiar velocity 
of the shell with radius Re centered on A3558. The plot of A'^ 
vs. Re depicts a step-like function where the value changes 
in discrete steps as observed clusters from the Table [T] be- 
come included in the subsample within Re - The discreteness 
of the limited sample is manifested in the sharp transitions 
exhibited in the radial profiles of the other quantities in 
Figure O The true profiles are smooth curves that roughly 
follow those calculated and are bounded by the ±1 — a error 
brackets approximately 68% of the time despite the sharp 
transitions. The uncertainty in the values of r, the distance 
between each cluster and A3558, shown in Table [1] is not 
included in the errors plotted in the figure. 

Since AMtot = 4-k R% pm{zssc) S, the plots of 5 and 
A Mtot vs. Re are analogous, however, they illustrate two 
different points. The radial profile of 5 shows a roughly 
power-law rise in the overdensity as we move toward the 
center of the SSC. The supercluster becomes quite nonlin- 
ear in the region interior to Re ~ 40 Mpc. The radial profile 
of AMtot, on the other hand, demonstrates that the mass 
only rises significantly through sharp jumps that correspond 
to the adding of additional clusters to the subsample. This 
can also be seen in the plot of Vp vs. Re where the profile 
between jumps is extremely fiat. Moreover, while the contri- 
bution to the peculiar velocity of the LG rises most rapidly 
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between Re ~ 30Mpc and Re ~ 45Mpc, Vp seems to have 
leveled off by Re ~ 50Mpc. 

The evolution of Vrad,p with Re is very mild with its 
mean varying by only a factor of ~ 1.5 over the full range of 
Re shown. In addition, the 1 — a errors on Vrad,p are large, 
comparable to the range over which Vrad,p varies. However, it 
is clear that the general trend is toward increasingly negative 
values of Vrad,p as we move toward the center of the SSC. 
This is what is expected from the spherical collapse model. 
In addition, by tracing the future of regions at different radii 
with the spherical collapse model (see ijS.ip . we determine 
that the SSC is only a closed system at Re ~ 10.3 given the 
mean of the density probability distribution and at 7?_b ~ 13 
if 1 — a fluctuations are considered. The mass enclosed within 
this radius is comparable to that of a rich X-ray cluster of 
~ IO^^Mq. The trajectories of spheres with Re ~ lOMpc, 
10.3 Mpc, 20Mpc, and 51 Mpc are plotted in Figure The 
mean values of the overdensity for each of these radii have 
been adopted in plotting the trajectories. 
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8 ROBUSTNESS OF RESULTS 

In ^ we showed how the matter overdensity in a region can 
be determined by fitting the mass function in an overdense 
region to the cluster mass function in the region. However, 
because of our limited sample, we only fit the normaliza- 
tion o f the mass function and relied on the lBarkana fc Loebl 
l|2004h hybrid model to correctly fit the slope. In this sec- 
tion, we investigate the robustness of our results and their 
dependence on our sample. In ^8.11 we consider the choice 
of obtained from the sample, while in i)8.2l we turn our 
attention to the slope of the mass function and ask whether 
different mass bins within the sample give consistent results. 



8.1 Dependence on A'^ 

The number of collapsed halos, their minimum mass, and 
their extent are the three key parameters that we extract 
from our cluster sample for use in our structure formation 
analysis. Though it is probably correct to identify the num- 
ber of X-ray luminous clusters with the number of collapsed 
halos, the two may not be exactly equivalent. In addition, 
there may be undiscovered objects obscured by galactic ab- 
sorption. As argued in SjS] the error in the mass estimates of 
each halo from the cluster X-ray properties can be expressed 
as an additional but small error in the number of halos above 
a fixed threshold. The same can be said for the extent of the 
SSC, denoted Re- Small errors in the physical position of a 
cluster with respect to the center of the SSC can be repre- 
sented as clusters being falsely included or excluded in our 
fixed selection of Re- We thus explore the dependence of 
our results on the assumed number of collapsed halos. 

Figure |4] shows the mean values of 5, b, M, Vrad,p, and 
Vp from Figure (TJ as well as the 1 — a errors for difi'erent 
choices of N - As A'' changes, Mmin and Re are held fixed 
at the values assumed throughout this paper. The results 
clearly do not change much as A'^ is varied over a relatively 
wide range of values. 



Figure 4. The dependence of the distributions for the nonlinear 
overdensity, 5, the total mass of the SSC, M, the cluster bias, b, 
the contribution of the SSC to the peculiar velocity of the LG, 
Vp, and the radial peculiar velocity of the outer edge of the SSC, 
^rad,py On the value of N. The square points represent the mean 
values of the probability distributions, while the error bars show 
the ±1 — (T spread. The vertical line denotes TV = 21. Mmin = 
1.75 X 10^** Mq and Re = 51 Mpc are held fixed for all values of 
N. 



8.2 Slope of the Mass Function 

Checking the fit of the mass function explicitly is equivalent 
to asking whether samples with higher mass thresholds give 
consistent results for the overdensity in the SSC. Since there 
are many more lower mass objects than higher mass ones, 
significant deviations in the number of higher mass objects 
may not affect the results obtained using all of the clusters. 
However, consistency between multiple mass bins shows that 
the results do not depend on the particular choices we made 
for our sample. 

The probability distributions for the matter overden- 
sity and the cluster bias are shown in Figure [5] calculated 
using each of four different mass thresholds. These distri- 
butions are not completely independent because the lower 
mass thresholds include all of the objects at higher mass. 
However, since there is a steep drop-off in the number of ob- 
jects with mass, the distributions are dominated by objects 
near the mass threshold. All four mass thresholds give con- 
sistent results for the matter overdensity in the SSC, with 
the distributions clustered around 5 ~ 0.8, slightly higher 
than the mean of the Mmin = 1-75 x 10^* Mq distribution. 
It is important to note that the mass threshold does not 
correlate with mean overdensity. This is just what we would 
expect if there is no bias in the number of objects with mass 
since, according to the model, a single value of the average 
linear matter overdensity should describe the entire mass 
function. 

On the other hand, we do expect a correlation between 
the bias and the mass of the objects in the sample. The figure 
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Figure 5. The probability distributions of the matter overden- 
sity, (S, (upper panel) and the cluster bias, b, (lower panel) ob- 
tained with four different cluster subsamples from Table [T] The 
solid curve shows the results from the entire sample of 21 clusters 
with masses above 1.75 X 10^* M0 shown in Figure [T] The dot- 
ted curve shows the results obtained using the 17 clusters with 
masses above 3.5 X lO^** Mq, while the 10 clusters with masses 
above 5.5 x 10^* Mq produced the results denoted by the dashed 
curve. Finally, the results denoted by the long-dashed curve were 
derived with only the 4 clusters with host halo masses above 
7.0 X IQi^Mo. The size of the region was fixed at = 51 Mpc 
for each distribution. 



clearly show this correlation as the mean of each bias dis- 
tribution increases with the mass threshold. This is a man- 
ifestation of the fact that more massive clusters are more 
strongly clustered than less massive ones. 



a radius of 51 Mpc, despite the significant mass in the region. 
The repulsive effect of the cosmological constant provides 
the extra push against gravity that will keep the region from 
ever collapsing. The outer shell is moving radially away from 
the center with a velocity only (437 ±87) kms~^ slower than 
the « 3600 km s~^ Hubble velocity at that radius. This pre- 
diction could be tested with better TuUy-Fisher data from 
SSC galaxies, or future CMB observations. 

We also calculated the contribution of the SSC to 
the peculiar velocity of the LG, Vp — (55 ± 13)kms^^. 
This value amounts to only (9.0 ± 2.1)% of the LG pe- 
culiar motion, muc h smal ler than a recent estimate by 
Kocevski fc Ebelind ( 2006 ) and more in agre ement with 
Ravchaudhurv et al.l ( 199ll ) , iReisenegger et all ((2000. ) , and 
Loeb fc NaravanI (|2007l V 

Despite the uncertainty in investigating the interior of 
the SSC using smaller subsamples of cluster, exploring the 
radial profiles of 5, Vp, and Vrad,p indicated convergence, i.e. 
that we have included all of the mass in the region in ex- 
cess of the universal average that contributes to the peculiar 
velocity of the LG. Moreover, while the 50 Mpc region as a 
whole is only mildly nonlinear, the interior of the SSC be- 
comes highly nonlinear. The excess mass becomes enough to 
bind the sphere with radius only ~ 10 Mpc. 

Finally, we showed that our results are robust to er- 
rors in our input parameters and the mass threshold of our 
cluster sample. Subsamples of clusters with different mass 
thresholds give consistent results for the overdensity in the 
region and demonstrate an expected increase in bias with 
mass. 
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9 CONCLUSIONS 

We have used the enhanced abundance of X-ray clusters to 
calculate the mass in th e SSC, based on the ePS-ST model 
IIBarkana fc Loebll2004 ). We constructed a sample of 21 X- 
ray luminous clusters within a radius of 51 Mpc in halos 
with masses above Mmin = 1.75 x 10^^ A/©. We then cal- 
culated probability distributions for the matter density con- 
trast of the SSC region, the cluster bias for our sample, 
and the mass of the SSC. We demonstrated that even mild 
values of the overdensity in a region can result in a signif- 
icant over-abundance of massive clusters. We found a mass 
of M = (4.4 ± 0.44) X 10^'^ Mq for the SS C. Our resuhs are 
in good agreement with previous results ( Reisenegger et al] 



I2OO0I : Ide Fihppis et al.ll2005l : iProust et al" 



20061 ). though we 



found that the cluster bias is probably higher than that es- 
timated bv lKocevski fc Ebelind (|2006l ). 

We then used the spherical collapse model to investigate 
the dynamics of the SSC. We found that the comoving size 
of the region has only collapsed by about 20% from its initial 
value. Moreover, we concluded that the SSC is not bound at 
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